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The abundance of primordial black holes is currently significantly constrained in a wide range of 
masses. The weakest limits are established for the small mass objects, where the small intensity of 
the associated physical phenomenon provides a challenge for current experiments. We used gamma- 
ray bursts with known redshifts detected by the Fermi Gamma-ray Burst Monitor (GBM) to search 
for the femtolensing effects caused by compact objects. The lack of femtolensing detection in the 
GBM data provides new evidence that primordial black holes in the mass range 5 x 10^^ - 10^° g 
do not constitute a major fraction of dark matter. 



7 I. INTRODUCTION 

8 Dark matter is one of the main and most challenging 

9 open problems in cosmology or particle physics, and a 

10 number of candidates for particle dark matter have been 

11 proposed over the years [ ] . An alternative idea that the 

12 missing matter consists of compact astrophysical objects 

13 was first proposed in the 1970s [ - ]. An example of 

14 such compact objects are primordial black holes (PBHs) 

15 created in the very early Universe from matter density 

16 perturbations. 

17 Recent advancements in experimental astrophysics, es- 

18 pecially the launch of the FERMI satellite with its un- 

19 precedented sensitivity has revived the interest in PBH 

20 physics [5, 6]. One of the most promising way to search 

21 for PBHs is to look for lensing effects caused by these 

22 compact objects. Since the Schwarzschild radius of PBH 

23 is comparable to the photon wavelength, the wave nature 

24 of electromagnetic radiation has to be taken into account. 

25 In such a case, lensing caused by PBHs introduces an in- 

26 terferometry pattern in the energy spectrum of the lensed 

27 object [ ]. The effect is called 'femtolensing' [ ] due to 

28 very small angular distance between the lensed images. 

29 The phenomenon has been a matter of extensive studies 

30 in the past [ ] , but the research was almost entirely the- 

31 oretical since no case of femtolensing has been detected 

32 as yet. Gould [ ] first suggested that the femtolensing 

33 of gamma-ray bursts (GRBs) at cosmological distances 

34 could be used to search for dark matter objects in the 

35 mass range 10^ — 10^^ g. Femtolensing could also be a 

36 signature of another dark matter candidate: clustered 

37 axions [ ]. 

38 In this paper, we present the results of a femtolensing 

39 search performed on the spectra of GRBs with known 

40 redshifts detected by the Gamma-ray Burst Monitor 

41 (GBM) on board the FERMI satellite. The non ob- 

42 servation of femtolensing on these bursts provides new 

43 constraints on the PBHs fraction in the mass range 



44 5 X 10-^^ — 10^^ g. We describe in details the optical 

45 depth derivation based on simulations applied to each 

46 burst individually. The sensitivity of the GBM to the 

47 femtolensing detection is also calculated. 

48 The paper is organized as follows: in Sec. II the basic 

49 equations for femtolensing and the calculation of lensing 

50 probability are given. Section HI describes the data sam- 

51 pie and simulations. In Sec. IV the results are presented, 

52 while Sec. V is devoted to discussion and conclusions. 



53 II. FEMTOLENSING 

54 A. Magnification and spectral pattern 

55 Consider a lensing event of a GRB by a compact ob- 

56 ject. The angular diameter distances from the observer 

57 to the lens, from the lens to the GRB source and from 

58 the observer to the source are Dql^ Dls^ and DqSi re- 

59 spectively. Coordinates are taken in the lens plane. The 

60 lens, with mass M, is located at the origin. The source 

61 position projected onto the lens plane is given by r^, the 

62 distance between the lens and true source position. The 

63 Einstein radius r^; is given by 

^2 ^ 4GM DqlDls 
c2 Dol + Dls 

64 where 



65 The image positions are given as usual by 

r± = l{rs±^Jrl + 4r%). (3) 
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66 The time delay St between the two images is given by 

cSt = V{ry,rs)-V{r^',rs), (4) 



2 



67 where V{r]rs) is the Fermat potential at the position r 

68 in the lens plane. One finds: 



cSt = 



2 \D 



LS 



DoL 



^^ln(--). (5) 



69 The phase shift between the two images is = E St/h^ 

70 where E is the energy of the photon. 

71 In the case of a point source, the amplitude contributed 

72 by the r± images is 



A± (X 



1 



(6) 



The magnification is obtained by summing the am- 
plitudes (6) and squaring, which gives 



\A[ 



1 



2cos(A0) 



1 



(7) 



The energy dependent magnification produces fringes 
in the energy spectrum of the lensed object. 

In principle, the finite size of the source and the relative 
motion of the observer, the lens and the source have to 
be taken into account. If the GRB is observed at a time 
texpi after the beginning of the burst, its size projected 
onto the lens plane is 



96 B. Lensing probability 

97 The lensing probability of gamma ray burst events is 

98 calculated in two steps. First, the optical depth r for 

99 lensing by compact objects is calculated according to 

100 the formalism of Fukugita et al. [ ]. The cosmologi- 

101 cal parameters used in the calculation are: a mean mass 

102 density Qm = 0.3 and a normalized cosmological con- 

103 stant = 0.7. The calculations are made for both the 

104 Friedmann-Lematre-Robertson- Walker (FLRW) and the 

105 Dyer-Roeder [12] cosmology. In our sample, the GRB 

106 redshift Zg is known. The lens redshift is assumed to 

107 be given by the maximum of the dr / dzL{zs) distribu- 

108 tion (see e.g. Fig. 5 of [11]). When r <C 1, the lensing 

109 probability p is given by p = r<j where a is the "lensing 
no cross-section" (see Chap. 11 of [ ]). 

111 In this paper, the cross-section is defined in the follow- 

112 ing way. Fringes are searched in the spectra of GRBs. 

113 These fringes are detectable only for certain positions rs 

114 of the source. The exact criteria for detectability will be 

115 given in section III C. The maximum and minimum po- 

116 sition of rs-, in units of ve are noted rs^min and rs^max- 
111 They are found by simulation and depend on the GRB 

118 redshift and luminosity. A minimum value of rs occurs 

119 because the period of the spectral fringes becomes larger 

120 than the GBM energy range at small rs- 

121 The femtolensing "cross-section" is simply 



sgrb 



DoL ct 



expl 



Dos r 

(c X 0.005 s) 



100 / V 2 



(8) 



82 where F is the Lorentz factor of the burst. Note that the 

83 Lorentz factor of GRBs is estimated to be in excess of 

84 100, so that sgrb given in Eq. (8) is overestimated. 

85 The ratio of sgrb to te is therefore 



Sgrb 



( ^expl \ 

' / M 
5GpcJ VlQiSg 



0.05 (^- 1)" 



Dos 



r 

100 



(9) 



Equation (9) shows that the finite size of the GRB can 
be in general safely neglected if texpi < 10 s. 

The Einstein radius crossing time tE is given by 



a = r 



S,max 



SjTnin 



(11) 



122 The lensing probability does not depend on the indi- 

123 vidual masses of lenses, but only on the density of com- 

124 pact objects ftco- In the optical depth calculation, an 

125 increase in the mass of the lenses is compensated by a 

126 decrease in the number of scatterers. Therefore, the con- 

127 straints for a given mass depend only on the cross section 

128 (J. 



129 III. DATA ANALYSIS 

130 In our analysis, we use a sample of GRBs with known 

131 redshifts. The selection of these bursts is described in 

132 Sec. Ill A. Then each burst is fitted to a standard spectral 

133 model, as explained in Sec. IIIB. Finally, the sensitivity 

134 of each burst to femtolensing is studied with simulated 

135 data. The simulation is described in Sec. IIIC. 



100 s 



rE 



c X 0.1 s 



300km/s 



(10) 



where v is the projected velocity of the source in the lens 
plane. Equation (10) shows that Ie ^ t^xpi under rea- 
sonable assumptions on the velocities. If so, the motion 
of the source in the lens plane can also be neglected. In 
the analysis of GRB spectra, it is thus assumed that the 
point source - point lens assumption is valid and that 
the source stays at a fixed position in the lens plane. 



136 A. Data selection 

137 The Gamma- Ray Burst detector (GBM)[ ] on-board 

138 the Fermi satellite consists of 12 Nal and 2 BGO scintil- 

139 lators which cover the energy range from 8 keV up to 40 

140 MeV in 128 energy bins. These detectors perform a whole 

141 sky monitoring. In the first two years of operation, GBM 

142 triggered on roughly 500 GRBs. In this paper, only the 

143 bursts with known redshifts have been investigated. The 
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144 initial sample consisted of 32 bursts taken from Gruber 

145 et al. [io] and 5 additional bursts from the GRB Coordi- 

146 nates Network (GCN) circulars [ ]. 

147 For 17 burst the amount of available data was not suf- 

148 ficient to obtain good quality spectra. The final sample 

149 thus consists of 20 bursts, which are listed in the Tab. I. 



195 1. The magnification (Eq. 7) as a function of the en- 

196 ergy is calculated for the given lens mass M and 

197 position of the source rs- 

198 2. This magnification is then convolved with the in- 

199 strumental resolution matrix to obtain magnifica- 

200 tion factors for each channel of the detector. 



150 B. Data processing and spectral analysis 

151 The GBM data are publicly available in the CSPEC 

152 format and were downloaded from the Fermi FSSC web- 

153 site [ ]. The CSPEC files contain the counts in 128 

154 energy channels binned in 1.024 s for all detectors. Only 

155 detectors with a minimal signal to noise ratio of 5.5 in 

156 each bin were selected for the analysis. 

157 Data were analyzed with the RMf it version 33pr7 pro- 

158 gram. The RMf it software package was originally devel- 

159 oped for the time-resolved analysis of BATSE GRB data 

160 but has been adapted to GBM and other instruments. 

161 For each detector with sufficient data, the background 

162 was subtracted and the counts spectrum of the first ten 

163 seconds of the burst (or less if the burst was shorter) was 

164 extracted. 

165 The energy spectrum was obtained with a standard 

166 forward-folding algorithm. Several GRB spectral mod- 

167 els such as a broken power law (BKN), Band's model 

168 (BAND) or a smoothly broken power law (SBKN) where 

169 considered. The femtolensing effect was added as a sepa- 

170 rate model. The magnification and the oscillating fridges 

171 where calculated according to Eq. (7), then multiplied 

172 with the BKN or BAND functions. 



173 C. Simulations 

174 The detectability of spectral fringes has been studied 

175 with simulated signals. The detectability depends on one 

176 side on the luminosity and the redshift of the bursts, and 

177 on the other side on the detector energy resolution and 

178 the data quality. The sensitivity of the GBM to the lens 

179 mass M depends strongly also on the energy range and 

180 resolution of the GBM detectors. When small masses 

181 are considered, the pattern of spectral fringes appears 

182 outside of the energy range. The large masses produce 

183 fringes with hardly detectable amplitudes and periods 

184 smaller than the energy bin size. 

185 Because the data quality and the background are 

186 not easily simulated, the detectability estimation is per- 

187 formed on real data. Namely, GRB events with known 

188 redshift are selected. Since the source redshift is known, 

189 the lens redshift is assumed to be the maximum value 

190 of dr /dzL{zs) as explained in Sec. II B. For a given ob- 

191 served GRB, the femtolensing signal depends thus only 

192 on 2 parameters: the lens mass M and the source posi- 

193 tion in the lens plane rs- The data are then processed as 

194 follows: 



201 3. The spectral signal is extracted from the data by 

202 subtracting the background. It is then multiplied 

203 by the corrected magnification. 

204 4. The background is added back. 

205 The detectability calculation can be illustrated with 

206 the luminous burst GRB090424592. The spectral data of 

207 this burst were first fitted with standard spectral models: 

208 BKN, SBKN and BAND. The GRB090424952 burst is 

209 best fitted with the BAND model. The fit has = 

210 78 for 67 degrees of freedom (d.o.f). The BAND model 

211 has 4 free parameters: the amplitude A, the low energy 

212 spectral index a, the high energy spectral index /3 and 

213 the peak energy Epeak V^]- 

214 The data are then modified by incorporating the spec- 

215 tral fringe patterns for a range of lens masses M and 

216 source positions rs- The simulated data and the corre- 

217 sponding femtolensing fit are presented in Fig. 1. Neither 

218 BKN nor BAND models are able to fit the simulated data 

219 (see Fig. 2). The values of rs are then changed until the 

220 of the fit obtained is not significantly different from 

221 the of the unmodified data. More precisely, the 

222 difference Ax^ should be distributed in the large sample 

223 limit as a x^ distribution with 2 degrees of freedom ac- 

224 cording to Wilk's theorem [ ]. The value Ax^ = 5.99, 

225 which corresponds to a x^ probability of 5% for 2 d.o.f, 

226 was taken as the cut value. The effect of changing rs on 

227 the femtolensing model is illustrated on Fig. 3 and 4. 

228 In Fig. 5 we show the maximum and minimum de- 

229 tectable rs for different lens masses. The maximum 

230 difference between rs^max and rs^min appears at M = 

231 5 X 10^^ g, which indicates the maximum of femtolensing 

232 cross-section. 



233 IV. RESULTS 

234 The 20 burst sample from Tab. I have been fitted with 

235 the standard BKN, BAND and SBKN models. The mod- 

236 els with the best x^ probability were selected and are 

237 shown on Tab. I. The bursts are well fitted by these stan- 

238 dard models, so that there is no evidence for femtolensing 

239 in the data. 

240 As explained in section II B, the lensing probability for 

241 each burst depends on the lens mass and on the rs^min 

242 and rs,max values. Since the sensitivity of GBM to fem- 

243 tolensing is maximal for lens masses of ~ 5 x 10^^ g 

244 (see Fig. 5), the values of rs,min and rs,max for each 

245 event were first determined at a mass M = 5 x 10^^ g 

246 by simulation. As explained in Sec. II B, the value of 
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FIG. 1. Simulated spectrum obtained with GRB 090424592. 
The spectrum was fitted with femtolensing+BAND model. 
The fit has — 92 for 73 d.o.f. The fit parameters are: 
A = 0.34di0.02phs-^ cm-2keV"\ E^eak = 173zbl2keV, a = 
-0.84 ±0.03 and /3 = -3.9 dz 7.5. The simulated femtolensing 
effect is caused by a lens at redshifts = 0.256 and a source 
at zs — 0.544. The simulated mass is M = 5 X 10^^ g and 
the mass reconstructed from the fit is 5.8 x 10^^ g. The source 
is simulated at position = 2. The position reconstructed 
from the fit is rs = 1-9. 
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FIG. 2. Simulated femtolensed spectrum fitted with the 
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parameters are: A = 0.5 ± 0.03 phs~^ cm~^ keV 
147zb5keV, a = -0.58±0.03 and /3 = -2.4d=0.1. 
model fit is almost indistinguishable from the BKN model fit 
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FIG. 3. The spectrum of GRB 090424592 using Nal detector 
n7, with the BAND and femtolensing fits superimposed. The 
parameters are rs = 1, 2, and lens mass 5 x 10^^ g. The 
models are convolved with the response matrix. 
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FIG. 4. The spectrum of GRB 090424592 using Nal detector 
n7. The BAND and femtolensing fits are superimposed. The 
parameters are rs — 3, 4, and lens mass 5 x 10^^ g. The 
excess at 33keV (K-edge) is an instrumental effect seen on 
many bright bursts. 



247 rs^min is Set by the period of the spectral fringes so that 

248 it is independent of the burst luminosity. The values of 

249 rs^max obtained are listed in Tab. 1. The lensing prob- 

250 ability is then calculated for both the FRLW and Dyer 

251 & Roeder cosmological models using each burst redshift, 

252 the most probable lens position and the values of rs^min 

253 and rs,max for the mass M = 5 x 10^^ g. The number 

254 of expected lensed bursts in the sample is the sum of the 

255 lensing probabilities. It depends linearly on flco- 

256 Since no femtolensing is observed, the number of ex- 



257 pected events should be less than 3 at 95% confidence 

258 level (C.L.). The constraints on the density of compact 

259 objects ftco is derived to be less than 4% at 95% C.L 

260 for both cosmological models. The values of the lensing 

261 probabilities for all the bursts in our sample assuming 

262 the constrained density of compact objects are shown in 

263 Tab. I. The limits at other lens masses are obtained by 

264 normalizing the Qco at M = 5 x lO^^g by the cross sec- 

265 tion a. The cross section is calculated using the Eq. 11 

266 and the values of rs^min and rs^max from Fig. 5. The 

267 limits on Vtco at 95% C.L. are plotted in Fig. 6. 
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271 dance of PBH is to constrain the density of compact ob- 

272 jects ^cO' Note that the hmits on the compact object 

273 abundance in the 10^^ — 10^^ g range obtained with mi- 

274 crolensing are at the 1% level. 

275 It is stated in Abramowicz et al. [ ] that the mass 
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276 range 10^^ g < Mbh < 10 g is virtually unconstrained. 

277 Under this range, the ^co is constrained by PBH evap- 

278 oration. Above this range, the constraints come from 

279 microlensing. The new idea by Griest et al. [ ] shows 

280 that the microlensing limit could be improved and get 

281 constraints down to 10^^ g with the Kepler satellite ob- 

282 servations. 

283 The FERMI satellite was launched three and a half 

284 years ago. Since then, almost 1000 of GRB were observed 

285 with the GBM detector. In many cases data quality is 

286 good enough to reconstruct time-resolved spectra. This 

287 unique feature is exploited in our femtolensing search by 

288 selecting the first few seconds of a burst in data analysis. 

289 Our limits were obtained by selecting only those bursts 

290 with known redshifts in the GBM data. This reduces the 

291 data sample from the 500 bursts detected in the first 2 

292 years to only 20. The constraints on ftco obtained at 

293 the 95% C.L. are shown on Fig. 6. These constraints 

294 improve the existing constraints by a factor of 4 in the 

295 mass range 5 x 10"^'^ - 10^^ g. 

296 After ten years of operation, the GBM detector should 

297 collect over 2500 bursts. Only a few of the bursts, say 

298 100, will have a measured redshift and sufficient spectral 

299 coverage. By applying the methods described in this pa- 

300 per, our limits will then improve by a factor of 5 reaching 

301 a sensitivity to density of compact objects down to the 

302 1% level. 



FIG. 6. Constraints on the fraction (or normalize density) 
of compact objects. The zones above the curves are excluded 
at the 95% confidence level. 
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TABLE I. The sample of 20 GBM GRBs used in the analysis. 













Fit to simulated data 










Name 




zs 


Fit to the data ^ 
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123/70 


103/68 


2 


0.346 


0.018 


0.020 


GRB 090424592 


0, 


.5440 


BAND 


78/67 


115/75 


97/73 


4 


0.256 


0.042 


0.046 


GRB 090510016 


0, 


.9030 


BKN 


62/66 


173/98 


158/96 


1.5 


0.406 


0.015 


0.016 


GRB 090618353 


0, 


.5400 


BAND 


59/72 


79.5/72 


66/70 


3 


0.254 


0.023 


0.025 


GRB 090926181 


2 


.1062 


BAND 


87/81 


105/81 


93/79 


4 


0.737 


0.413 


0.423 


GRB 091003191 


0, 


.8969 


BKN 


93/94 


105/94 


94/92 


3 


0.400 


0.059 


0.064 


GRB 091020900 


1, 


.7100 


BKN 


74/69 


116/69 


90/67 


2.5 


0.667 


0.119 


0.124 


GRB 091127976 


0, 


.4900 


BAND 


78/74 


84/74 


76/72 


4 


0.240 


0.034 


0.037 


GRB 091208410 


1, 


.0630 


BAND 


55/55 


115/59 


56/57 


2.5 


0.457 


0.055 


0.059 


GRB 100414097 


1, 


.3680 


BKN 


65/61 


97.5/68 


86/66 


2.5 


0.560 


0.084 


0.089 


GRB 100814160A 


1, 


.4400 


BKN 


42/40 


138/42 


128/40 


2 


0.590 


0.058 


0.061 


GRB 100816009 


0, 


.8049 


BKN 


73/52 


95/50 


66/48 


2.5 


0.360 


0.034 


0.037 


GRB 110731465 


2 


.8300 


SBKN 


72/64 


125/64 


97/62 


3 


0.877 


0.347 


0.344 



^ Fit has been performed using only the photons arrived in less than 10s from the beginning of the burst. 
^ for assumed ^co = 0.037 
^ for assumed flco = 0.041 



